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ABSTRACT 

A parameter study of the spectral evolution of a typical post-AGB star, with particular 
emphasis on the evolution of the infrared colours, is presented. The models are based 
on the latest evolutionary tracks for hydrogen burning post-AGB stars. For such tracks 
the evolutionary rate is very dependent on the assumed mass loss rate as a function 
of time. We investigate this effect by modifying the mass loss prescription. The newly 
calculated evolutionary rates and density distributions are used to model the spectral 
evolution of a post-AGB star with the photo-ionization code CLOUDY, including dust in 
the radiative transfer. Different assumptions for the dust properties and dust formation 
are considered. It is shown that by varying these parameters in a reasonable way, 
entirely different paths are followed in the IRAS colour-colour diagram. First of all, 
the effects of the evolution of the central star on the expanding dust shell can not be 
neglected. Also the dust properties and the definition of the end of the AGB phase 
have an important effect. The model tracks show that objects occupying the same 
location in the IRAS colour-colour diagram can have a different evolutionary past, 
and therefore the position in the IRAS colour-colour diagram alone can not a priori 
give a unique determination of the evolutionary status of an object. An alternative 
colour-colour diagram, the K-[12] vs. [l2]-[25] diagram, is presented. The tracks in this 
diagram seem less affected by particulars of the grain emission. This diagram may be 
a valuable additional tool for studying post-AGB evolution. 
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1 INTRODUCTION 

The transition phase between the Asymptotic Giant Branch 
(AGB) and planetary nebulae (PNe) has gained much atten- 
tion over the last decade. AGB stars lose mass fast and get 
obscured by their circumstellar dust. As the star leaves the 
AGB, its mass loss rate decreases significantly and the star 
may become suSiciently hot to ionize its circumstellar mate- 
rial and be observable as a PN. During the transition from 
the AGB to the PN phase (the post-AGB or proto-planetary 
nebula phase) the dust shell created during the AGB moves 
away from the central star and becomes optically thin after 
a few hundred years; the obscured star becomes observable. 
The transition time from the AGB to the PN phase is esti- 
mated to be a few thousand years (e.g. Pottasch 1984). 

Whereas PN are relatively easy to find because of their 
rich optical emission line spectra, post-AGB stars have more 
inconspicuous spectra and are therefore much harder to find. 
The number of post-AGB stars only started to become large 



after the IRAS mission, that was successful in detecting ob- 
jects surrounded by circumstellar dust. 

Several samples of post-AGB stars are presented in the 
literature (e.g. Volk & Kwok 1989, Hrivnak, Kwok & Volk 
1989, van der Veen, Habing & Geballe 1989, Trams et al. 

1991, Oudmaijer et al. 1992, Slijkhuis 1992). Most of these 
objects are stars with supergiant-type spectra, surrounded 
by dust shells. These samples of post-AGB stars have in 
common that the original criteria which were employed to 
find them implicitly made assumptions on the spectral en- 
ergy distribution (SED) of post-AGB stars. Some authors 
used criteria on the IRAS colours, because post-AGB stars 
were expected to be located in a region in the IRAS colour- 
colour diagram between AGB stars and PN (e.g. Volk & 
Kwok 1989, Hrivnak et al. 1989, van der Veen et al. 1989, 
Slijkhuis 1992). Other authors loosened this criterion and 
searched for objects in the entire colour-colour diagram, but 
with an additional criterion that the central star should be 
optically visible (e.g. Trams et al. 1991, Oudmaijer et al. 

1992, Oudmaijer 1996). 
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Such samples are subject to selection effects, so the ob- 
jects that have been selected do not have to be representa- 
tive for the entire population of post-AGB stars. In order 
to understand these selection effects and to obtain a handle 
on the kinds of objects that could have been missed, it is 
useful to investigate the spectral energy distribution from 
a theoretical point of view by following the spectral evolu- 
tion of a post-AGB star with an expanding circumstellar 
shell. Moreover, this type of study allows one to investigate 
and understand the processes that occur in the circumstellar 
shell during the transition. 

Several such studies have been published. Most of these 
studies focus on the expanding dust shell with a dust radia- 
tive transfer model. Authors like e.g. Siebenmorgen, Zijl- 
stra & Kriigel (1994), Szczerba & Marten (1993), Loup 
(1991), Slijkhuis & Groenewegen (1992) and Volk & Kwok 
(1989) performed calculations describing the evolution of 
the circumstellar dust shell. Work concerning the evolution 
of a star with an expanding shell has also been performed 
with photo-ionization codes (Volk 1992), and hydrodynam- 
ical models (Frank et al. 1993, Mellema 1993, Marten & 
Schonberner 1991). None of these models include dust, ex- 
cept for Volk (1992) who used the output of the photo- 
ionization code CLOUDY (Ferland 1993) as input for a dust 
model. 

Our objective is to investigate the spectral evolution of 
a hydrogen burning post-AGB star with a photo-ionization 
model containing a dust code. The aim of this work is 
twofold: 

Firstly we investigate the processes in the circumstellar en- 
velope. The emphasis in this paper will be on the infrared 
properties of this shell. For this, both the expansion of the 
shell and the evolution of the central star have to be taken 
into account. 

Secondly, we investigate the influence of certain assumptions 
on the evolutionary timescales. Only few post-AGB evolu- 
tionary grids have been published, never giving a fine grid for 
the coolest part of the evolution. The original Schonberner 
tracks (1979, 1983) presented only a limited number of time 
points of the evolution, while Vassiliadis & Wood (1994) 
omit the phase between 5000 K and 10 000 K altogether. 
The predicted timescales that are available are calculated 
with a pre-defined end of the AGB and assumed post- 
AGB mass loss rates. These choices can influence the post- 
AGB evolutionary timescales considerably, as was already 
demonstrated by Trams et al. (1989) and Gorny, Tylenda 
& Szczerba (1994). The situation has changed now with the 
results of Blocker (1995a, b), who calculated new evolution- 
ary sequences, and made extensive tables available describ- 
ing certain key parameters during the post-AGB phase. The 
published relation between the envelope mass and effective 
temperature allows one to construct detailed timescales us- 
ing one's own mass loss prescriptions during the (post-) AGB 
evolution, which makes the model results less dependent on 
the mass loss formulation that was used by Schonberner 
(1979, 1981, 1983). 

In this paper first we describe the method to calculate 
evolutionary timescales and the adopted mass loss prescrip- 
tions. The results of Blocker (1995a, b) are used as a basis 
to create synthetic evolutionary tracks, and some aspects 
of the evolutionary timescales that are predicted are dis- 
cussed. Next we describe the photo-ionization code CLOUDY 



that was used and the assumptions that were made to con- 
duct the study of the spectral evolution of post-AGB stars. 
We then present the first results of a parameter study of a 
typical post-AGB object, based on the 0.605 M© track from 
Blocker (1995b). 



2 THE CENTRAL STAR EVOLUTION 

Many stellar evolutionary models are presented in the recent 
literature, but only two groups calculate the AGB quanti- 
tatively including mass loss. These are Vassiliadis & Wood 
(1993, 1994) and Blocker & Schonberner (1991) and Blocker 
(1995a,b). Both groups calculate the evolution of a star from 
the main sequence through the red giant phase to the white 
dwarf stage. There are slight differences in the core mass - 
luminosity relations and the use of a different initial - fi- 
nal mass relation, but their results show qualitatively the 
same behaviour in the evolution of stars in the HR diagram. 
The main differences between the models are the mass loss 
prescriptions on the AGB. 

The AGB mass loss rates in the formulation of Vassil- 
iadis & Wood are derived from the mass loss - pulsation 
period (M-P) relation given by Wood (1990). For periods 
longer than 500 d for low mass stars, and larger pulsation pe- 
riods for higher mass objects, a maximum value of the AGB 
mass loss rate is invoked (of the order of 10^^ Meyr^^). 
Blocker used a mass loss rate which is dependent on the lu- 
minosity of the star. He fitted the results of Bowen's (1988) 
theoretical study of mass loss in Mira variables. Basically 
this is the Reimers mass loss (Reimers 1975) multiplied by 
a luminosity dependent factor. These mass loss rates are 
not limited by a maximum value. Since the adopted AGB 
mass loss rates by the above authors differ, some striking 
differences between the results of their calculations exist. 
The M-P relation of Wood (1990) is questioned by Groe- 
newegen & de Jong (1994). Using the synthetic evolutionary 
model of Groenewegen & de Jong (1993), they were able to 
fit the luminosity function of carbon stars in the LMC with 
the Bowen mass loss adopted by Blocker (1995a), but not 
with the Vassiliadis & Wood mass loss rates. 

The mass loss rates do not only govern the evolution 
on the AGB. During the post-AGB phase, the mass loss 
rates also have a drastic infiuence on the timescales of the 
AGB - PN transition. Since the temperature of the star for 
a given core mass is determined by the mass of the stellar 
envelope, larger mass loss rates will cause the star to evolve 
to higher temperatures more quickly and can therefore de- 
crease the timescale of the transition strongly. Trams et al. 
(1989) showed that when the adopted post-AGB mass loss 
rates of a 0.546 M© star are raised by a factor of 5 to 10 to a 
value of 10~^ M© yr~^, the transition time from the AGB to 
the PN phase is shortened from 100 000 yr to only 5000 yr. 
This would make a low mass star readily observable as a PN, 
while the (longer) predicted timescales prevent such objects 
to become an observable PN since the circumstellar shell 
would have moved far away from the star and would have 
dispersed into the interstellar medium long before the star 
emits a sufficient amount of ionizing photons. In order to 
study the evolutionary timescales of post-AGB stars, it is 
important to understand post-AGB mass loss better. 
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2.1 Determination of the evolutionary timescales 

The evolutionary timescales of hydrogen burning post-AGB 
stars are determined by the core mass and the mass loss rates 
of the star. In this section we present the computational de- 
tails of this procedure, compute the evolutionary timescales 
and investigate several aspects of these timescales. We use 
the tracks by Blocker (1995a,b) to determine the evolution- 
ary timescales. Blocker calculated several tracks for his in- 
vestigation of the evolution of stars on the AGB and beyond. 
The tracks include a complete calculation of a 3 M0 object 
that ends as a 0.605 M© white dwarf and a 4 M© (0.696 M©) 
sequence. To these sequences, the 1.0 Mq (0.565 Mq) track 
of Schonberner (1983) was added. We start with a recapit- 
ulation of the mass loss laws we have used (see also Blocker 
1995b). 



2.1.1 The mass loss laws 

During the red giant branch, the mass loss is parameterized 
by the Reimers mass loss rate Mr given in equation (|l|), 
where the scaling parameter 77 is set to 1 for all tracks. 



M./(Meyr-)^4xlO-,i^l/^|)^ 



(1) 



With L, R and M the luminosity, radius and mass of the 
star. Since AGB stars lose mass at a faster rate than red 
giants, a different formulation is adopted for the AGB mass 
loss rates. To this end. Blocker fitted the numerical results 
of Bowen, as given in equation (Gl). 

MBi/(M0yr-^)= (2) 

4.83 X 10-8(Mzams/Mo)-2-1(L/Lo)'''Mr 

MB2/(Moyr-i) =4.83 x 10"' (M/Mq)-''-'(L/Lq)"'"Mr 

The difference between AIbi and Mb2 is the division by the 
current mass of the star instead of the initial mass of the 
star. The timescales for the 0.565 M© track of Schonberner 
were recalculated by us using the Monv-reff relation given by 
Schonberner but with the mass loss prescriptions of Blocker. 
For all tracks the Mbi prescription was used, except for the 
0.605 M0 track where Mb2 was used, which results in larger 
AGB mass loss rates. 

The mass loss rates can reach values of the order of 
10"* M0yr"^ to 10"^ Mgyr'^ and it is clear that if this 
mass loss would last for a long time, the star would evapo- 
rate. Thus, an end to the AGB (-wind) has to be invoked. 
Blocker used the pulsation period to define the end of the 
high mass loss phase; he assumed that AGB stars pulsate in 
the fundamental mode, where the period can be calculated 
using equation (H) (Ostlie & Cox 1986). 

lg(Po/d) = -1.92 - 0.73 lg(M/M0) + 1.86 lg(i?/R0) (3) 

When the central star has reached an inferred pulsation pe- 
riod of Pa = 100 d (which occurs when the star has a sur- 
face temperature somewhere between roughly 4500 K and 
6000 K) the Bowen mass loss stops. When the star has sub- 
sequently reached an inferred pulsation period of Pb = 50 d 
the post-AGB mass loss starts; in between the mass loss 
rates are connected by a smooth transition. Hence the fol- 
lowing definitions will be used in the remainder of the paper: 



the AGB phase is that part of the evolution where the pul- 
sation period is greater than Pb, the (AGB) transition phase 
is that part of the evolution where the pulsation period is in 
between Pa and Pb, and finally the post-AGB phase is that 
part of the evolution where the pulsation period is smaller 
than Pb. 

One should realize that post-AGB stars with pulsation 
periods larger than 50 d exist (e.g. HD 52961 with a period of 
72 d; Waelkens et al. 1991, Fernie 1995). Thus the pulsation 
period recipe used by Blocker should only be regarded as an 
approximate parameterization of the end of the AGB. 

For lower temperatures the post-AGB mass loss is given 
by the Reimers law (llf). Since the Reimers mass loss is 
proportional to TcS~^ for a constant luminosity and stel- 
lar mass, the post-AGB mass loss rates decrease during the 
evolution of the object. A radiation driven wind will take 
over when the star has reached temperatures above approx- 
imately 20 000 K. The mass loss rate for this wind, based 
on Pauldrach et al. (1988), is given by equation M). Hence 
at any stage of the post-AGB evolution either equation (uh 
or (W) is used, whichever of the two yields the biggest mass 
loss. 



McPN/(M0yr-') = 1.29 x 10^'" (L/Lq)' 



(4) 



During the evolution on the AGB and beyond, the en- 
velope mass is reduced due to two processes. In the outer 
parts of the envelope, mass is lost through a wind, and at 
the bottom of the envelope, mass is diminished by hydrogen 
burning. From Trams et al. (1989) we find the mass loss due 
to hydrogen burning; 



MH/(MQyr"^) = 1.012 x 10~" (L/Lq)A:; 



(5) 



With Xe the hydrogen mass fraction in the envelope (70 %) 



2.1.2 The evolutionary timescales 

The evolutionary timescales for hydrogen burning post-AGB 
stars depend firstly on the core mass and secondly on the 
mass loss rates. The evolutionary timescales can be calcu- 
lated relatively easy by making use of the fact that there 
exists a unique relation between the envelope mass and the 
stellar temperature for every core mass. One can calculate 
the evolutionary timescales by combining this relation with 
the mass loss prescriptions given in the previous section. 
Dr. Blocker kindly provided us with tables from which the 
Manv-Tcs relation could be reproduced. 

The timescale for the envelope depletion is given by: 



AAfe, 



= A/wind + Mh 



(6) 



At 

Using the Monv-roff relation this expression can be trans- 
formed into an expression for the evolutionary rate ATcg/At 
of the central star 



At 



AT, 



AAL, 



cff Ai\4, 

X 



At 



Mwind + Mh 
di\'ic„v 



(7) 



Integrating this expression yields the actual timescales for 
the evolution of the central star. 

In Fig. |l|some relations are presented for the 0.565 M0, 
0.605 M© and 0.696 M© tracks. The upper panel shows the 
envelope mass as function of photospheric temperature. The 
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Figure 1. The Menv— Teff relations, mass loss rates and result- 
ing evolutionary rates for the 0.565 M-q (dotted line), 0.605 M© 
(solid line) and 0.696 Mq (dashed line) tracks. The envelope mass 
is given in solar masses, the mass loss rate in solar masses per 
year and the evolutionary rate in kelvin per year. Note that the 
0.605 Mq track has a larger AGB mass loss rate. For this partic- 
ular track, the B2 mass loss rate was used. 



second panel presents the mass loss rates as function of tem- 
perature, and the third panel shows the evolutionary rate 
in kelvin per year as function of temperature. In the sec- 
ond panel it is visible that the transition phase occurs at 
higher temperatures for larger core masses with the recipes 
described above. 

The evolutionary rate, mass loss rate and envelope mass 
are related to each other in the following way. The evolution- 
ary rate in terms of increase in temperature per year is slow 
on the steep part of the Mcnv-Teff relation, and more rapid 
on the shallow part. Larger mass loss rates imply of course 
more rapid changes in the stellar temperature. These effects 
are visible in Fig. hi: at low temperatures in the post- AGB 
phase the evolutionary rates are smallest, while for higher 
temperatures the evolutionary rate increases. The minima 
around lg(rcff/K) ~ 3.7 to 3.9 correspond to the onset of 
the smaller post-AGB mass loss, slowing down the evolution 
which accelerates later on the shallow part of the Menv-rcff 
relation. 

Interestingly, the minimum in the evolutionary rate oc- 
curs right after the end of the AGB transition phase. The 
net increase in effective temperature in kelvin per year is the 
slowest for all tracks just after the start of the post-AGB evo- 
lution, which is when the temperature of the objects corre- 
spond to G or F spectral types. The increase in temperature 
is less than 1 Kyr"^ for the 0.565 M© and 0.605 M© tracks. 



How does this evolutionary rate compare with the obser- 
vations? Fernie & Sasselov (1989) calculated the possible 
increase in temperature for UU Her stars. From the absence 
of a change in pulsation period of UU Her, 89 Her and HD 
161796, they place an upper limit of 0.5 Kyr^^ on the tem- 
perature increase. Their conclusion was that these objects 
can not be post-AGB stars because the evolutionary rates 
should be much higher. However, this may be the case when 
one averages over the entire post-AGB temperature span, 
but the observed lack of temperature increase is consistent 
with the predictions for this temperature range, as was al- 
ready shown by Schonberner & Blocker (1993). Therefore a 
post-AGB nature for the UU Her stars can not be excluded 
on this basis. 



2.2 The influence of mass loss on the evolutionary 
timescales 

It is evident from the above that the value of the post-AGB 
mass loss rate has an important effect on the evolutionary 
timescales. However, the mass loss rate is not known ob- 
servationally for cool post-AGB stars. The usual tracer of 
mass loss, Ha emission, which is often observed in the spec- 
tra of post-AGB objects is likely to be the result of stellar 
pulsations (see the discussion by Oudmaijer & Bakker 1994 
and Lebre et al. 1996). A possible tracer of mass loss in cool 
post-AGB stars is the CO first-overtone emission at 2.3 /im 
(Oudmaijer et al. 1995), but, given this is true, the mass loss 
rates still have to be determined. 

The lack of theoretical and observational values for 
post-AGB mass loss rates forced Schonberner and Blocker 
to resort to the heuristic Reimers law for the cool part of the 
post-AGB evolution. As an illustration of the effect of the 
post-AGB mass loss rates on the evolutionary timescales, we 
have calculated these timescales for three core masses with 
the post-AGB mass loss rate at 0, 1, 5 and 10 times the 
standard post-AGB value (indicated as OxpAGB etc.). The 
results for the 0.565 M©^ 0.605 M© and 0.696 Mq tracks 
are presented in Table hi, where the timescales since the 
end of the transition mass loss phase are given. The in- 
crease of the post-AGB mass loss rates indeed decreases the 
timescales of the evolution, confirming the results of Trams 
et al. (1989) and Gorny et al. (1994). The OxpAGB mass loss 
rate effectively determines the slowest possible evolutionary 
speed, since the only mass loss is through hydrogen burning. 
On average, the difference in speed between OxpAGB and 
IxpAGB is roughly a factor of two to three. 



2.3 Distribution over spectral type 

The availability of the evolutionary rates allows us to inves- 
tigate the predicted distribution over spectral type. Oud- 
maijer, Waters & Pottasch (1993) and Oudmaijer (1996) 
used the coarse grid of the Schonberner (1979, 1983) tracks 
and found that a star spends by far most of the time as a 
B-type star. Only for the 0.644 Mq track half of the time is 
spent in the G phase, and somewhat less as a B star, while 
almost no time is spent as an F or A star. 

One might not expect that an object would spend a 
large fraction of the time as a B star, because the evolution- 
ary rates are largest for B spectral type (Fig. fll). This can 
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Table 1 . Evolutionary timescales for different values of post-AGB 
mass loss rates; the post-AGB phase starts at Pb = 50 d. 



-'W'pAGB OxpAGB 
%S (K) 



IxpAGB SxpAGB lOxpAGB 



Af„ 



0.565 Mq, L = 3891 L^ 



5 426 














5 500 


312 


76 


18 


9 


6 000 


1113 


286 


70 


36 


7000 


2810 


815 


207 


107 


10 000 


3974 


1306 


349 


182 


15 000 


5172 


2092 


651 


353 


20 000 


6383 


3089 


1141 


653 


25 000 


7544 


4156 


1768 


1068 


30 000 


8709 


5292 


2526 


1604 




A/coro = 0.605 Mq,L-- 


= 6310 L© 




6 042 














6 500 


1671 


458 


113 


58 


7000 


2510 


715 


178 


92 


10 000 


3122 


948 


241 


125 


15 000 


3487 


1181 


325 


171 


20 000 


3822 


1458 


452 


248 


25 000 


4154 


1769 


625 


358 


30 000 


4481 


2093 


832 


501 




Afcorc = 0.696 Mq, L = 


11610 Lq 




6 846 














7000 


94 


29 


7 


4 


7500 


192 


62 


15 


8 


10 000 


245 


84 


22 


11 


15 000 


287 


111 


30 


16 


20 000 


309 


129 


38 


21 


25 000 


331 


151 


50 


28 


30 000 


354 


174 


65 


38 



be understood however when we consider the large range of 
temperatures that corresponds to spectral type B: roughly 
between 10 000 K and 30 000 K. In contrast, A stars only 
have a temperature range between approximately 7500 K 
and 10 000 K. The large evolutionary rates multiplied by 
the temperature range then result in a longer time spent as 
a B star during the post-AGB evolution. The large fraction 
of time that is spent as a G star in the 0.644 Mq track is ex- 
plained by the steep M^-nv-T^s relation for temperatures less 
than approximately 6000 K, as can be deduced from equa- 
tion (rTI). The distribution over spectral type for the tracks 
presented here with IxpAGB mass loss is calculated using 
the conversion from effective temperature to spectral type 
listed by Straizys & Kuriliene (1981). These are TefF(AOI) = 
9800 K, T;ff(FOI) = 7400 K and rcff(GOI) = 5700 K. 

To allow for a comparison with the distributions pre- 
sented by Oudmaijer et al. (1993) we will assume in this 
section that the post-AGB phase starts when Tefi = 5000 K 
and ends when TcS — 25 000 K. The distributions can be 
easily obtained by calculating the time spent as a G star 
(5000 K - GO), F star (GO - FO) etc., and subsequently di- 
viding these numbers by the total time spent as a post-AGB 
star. The resulting distributions are plotted in Fig. H. For 
comparison the distribution over spectral type of the sample 
of 21 post-AGB objects in the list of Oudmaijer et al. (1992) 



a. 




Figure 2. Distribution over spectral type. Upper left panel; the 
post-AGB candidate stars in the sample of Oudmaijer et al. 
(1992), the other panels show the predicted distribution for the 
0.565 Mq, 0.605 Mq and 0.696 Mq tracks. 



is given. The observed distribution peaks at F, while no B 
stars are found at alLn 

One should realize that planetary nebula central stars 
with temperatures below 25 000 K are observed: e.g. 
IRAS19336-0400 with Tch = 23 000 K (Van de Steene, Ja- 
coby & Pottasch 1996, Van de Steene & van Hoof 1995). 
Hence one could also assume an upper limit of 20 000 K for 
the post-AGB regime. However, the distribution over spec- 
tral type is not very different in this case, the fraction of 
B-type stars will be lower by approximately 10 %. 

The 0.605 Mq distribution is different from the results 
presented by Oudmaijer et al. (1993) for the 0.598 Mq 
Schonberner track. In the present plot, the distribution 
peaks at F, while the 0.598 Mq distribution peaks at G. 
The distribution for the same track is different in Oud- 
maijer (1996), since there the evolutionary timescales were 
shortened by 1000 yr 'in order to make the old 0.598 Mq 
track consistent with the new calculations' (Marten & 
Schonberner 1991). This resulted in a distribution that 
strongly peaks at spectral type B. 

The main differences between the old 0.598 Mq and new 
0.605 Mq calculations are the progenitor mass of the star 
(1 Mq and 3 Mq respectively), the mass loss prescriptions 
and the definition of the end of the AGB. In the Schonberner 
calculations, no transition wind was assumed: the AGB mass 
loss would abruptly change into a post-AGB wind at 5000 K. 
In the Blocker calculations, a transition wind is assumed be- 



* The few B type objects in the sample of Oudmaijer et al. 
(1992) appear to have low effective temperatures. HR 4049 and 
HD 44179 (the central star of the Red Rectangle) are listed in 
the literature as a B star, but abundance analyses showed that 
the effective temperature is lower than typical for a B star: about 
7500 K (late A or early F type; van Hoof et al. 1991, Waelkens 
et al. 1992). The reason that these objects were classified as B, 
instead of late A stars, is the extreme metal deficiency of these 
stars, so that few metallic lines are present in the spectra. Thus 
the observed spectra mimic the spectra of hot objects. 
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Table 2. The results of the Kolmogorov-Smirnov test comparing 
various evolutionary tracks with the observed post-AGB sample. 
The tests have been performed using a range of temperatures for 
the post-AGB regime of 5000 K to 20 000 K (second column) or 
alternatively 5000 K to 25 000 K (third column). 



track probability probability 
Mq 20 000 K 25 000 K 



0.565 
0.605 
0.696 



4.5xl0~* 
2.4x10-2 
6.7x10-2 



2.8x10-*^ 
6.4x10"" 
1.3x10-1 



tween 5000 K and 6000 K for the 0.605 Mq track. This im- 
plies s/iorter timescales for the G type phase of the 0.605 Mq 
track with respect to the older calculations, but does not af- 
fect the time spent as F, A or B stars, where for both tracks 
a Reimers post-AGB mass loss is assumed. The only expla- 
nation we can find for the apparent difference in evolution- 
ary speed is that the dependence of the temperature on the 
envelope mass of the 0.605 Mq track is different from the 
0.598 Mq track. Apparently, the slope of the Menv-Toff re- 
lation can differ significantly for stars with a different evolu- 
tionary past, even when the resulting core masses are nearly 
identical. 

The predicted distributions in Fig. ti show large differ- 
ences. For the smallest core mass, most of the time is spent 
as a B star. For the 0.605 Mq case the lifetime is distributed 
evenly over the F and B phase. The fastest track with a core 
mass of 0.696 Mq shows a minimum at A, while the rest of 
the lifetime is spread evenly over G, F and B. The 0.605 Mq 
distribution nicely reproduces the observed peak at F. How- 
ever, all tracks predict many more B stars than observed. 

In order to compare the observed and predicted 
distributions more quantitatively, we have performed a 
Kolmogorov-Smirnov test using both an upper limit of 
20 000 K and 25 000 K for the post-AGB regime. The re- 
sulting probabilities are given in Table 0. From these tests it 
appears that 0.565 Mq track is an unlikely model for the ob- 
served post-AGB sample. On the other hand, the 0.605 Mq 
and the 0.696 Mq tracks can not be excluded. Given the fact 
that stars on the 0.696 Mq track evolve much faster than 
on the 0.605 Mq track, they would constitute only a small 
fraction of the total number of observable post-AGB stars. 
We will adopt the 0.605 Mq track to describe the post-AGB 
evolution in the remainder of this paper. 

This exercise shows that the predictions of the distri- 
bution over spectral type are subject to large uncertainties, 
depending both on the A'lanv-Tag relation and the assumed 
mass loss prescription. However, regardless what assump- 
tions are made, one would expect a fair number of B-type 
post-AGB stars, which are not observed in the sample de- 
picted in Fig. pi Hence this discrepancy remains unresolved. 
Reversely it can be stated that when a larger sample of 
post-AGB stars with reliable temperature determinations 
becomes available, the procedure described here can be a 
very effective means for testing evolutionary tracks. 



3 THE MODEL 

In order to calculate the spectral evolution of post-AGB 
stars, we used the photo-ionization code CLOUDY version 
84.12a (Ferland 1993). Some modifications have been made 
to the code to facilitate the computations. The most impor- 
tant change was the introduction of several broadband pho- 
tometric filters, including the Johnson and the IRAS filters. 
The in-band fluxes for these filters were calculated by fold- 
ing both the spectral energy distribution and the emission 
line contribution with the filter passband. Internal extinc- 
tion due to continuum opacities was included both in the 
continuum and the line contribution. 

A dust model written by P.G. Martin was already in- 
cluded in the original code. For our modeling we used the 
grain species labelled 'ISM Silicate' and TSM Graphite'. The 
optical constants were taken from Martin & Rouleau (1991). 
The absorption and scattering cross sections were calculated 
assuming a standard ISM grain size distribution (Mathis, 
Rumpl & Nordsieck 1977). All calculations were done as- 
suming a dust-to-gas mass ratio of 1/150. For the chemical 
composition of the gas we assumed the abundances given in 
AUer & Czyzak (1983), supplemented with educated guesses 
for elements not listed therein (as given in cloudy). 

The original code only allowed for the computation of a 
model with a constant dust-to-gas ratio throughout the en- 
tire nebula. However, the density profiles we calculate extend 
from the stellar surface outward. It is not realistic to assume 
that dust is present near the stellar surface and therefore we 
introduced new code in CLOUDY to solve this problem. This 
enabled the dust to exist only outside a prescribed radius 
or, alternatively, only in those regions where the equilib- 
rium temperature of the dust would be below a prescribed 
sublimation temperature. These prescriptions work as a bi- 
nary switch: at a certain radius either no dust or the full 
amount is present. In those regions where dust exists, the 
dust-to-gas ratio is assumed to be constant. 

Two models for the dust formation are adopted. In the 
first model it is assumed that dust is only formed in the AGB 
wind, hence in material that was ejected before the stellar 
pulsation period reached Ph ■ We will call this the AGB-only 
dust formation model. In the second model it is assumed 
that the dust formation continues in the post-AGB wind. 
We will call this the post-AGB dust formation model. Due 
to limitations of the code, which we discuss below, we will 
only investigate the spectral evolution after the post-AGB 
phase has started. This implies that for the AGB-only dust 
formation model, the inner dust radius is already at a dis- 
tance from the central star and the equilibrium temperature 
of the grains is always below the sublimation temperature. 
In the post-AGB dust formation model we assume that dust 
only exists in those parts of the nebula where the equilibrium 
temperature of the grains is below the sublimation temper- 
ature. The assumed values for the sublimation temperature 
are 1500 K for graphite and 1000 K for silicates. 

It should be noted that we only try to model dust for- 
mation and not the destruction of grains by the stellar UV 
field or shocks. Especially in the AGB-only dust formation 
model the grains at the inner dust radius are always exposed 
and it is expected that they eventually will be destroyed. 
However, little is known about grain destruction, and the 
rates at which this destruction occurs are very uncertain. In 
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the case of continuing dust formation in the post-AGB wind 
this problem can be expected to be of lesser importance since 
there is a constant supply of new grains shielding the older 
grains. 

The radiative transport in CLOUDY is treated in one 
dimension only, i.e. the equations are solved radially out- 
wards. This assumption makes the code unsuitable to com- 
pute models of nebulae with significant amounts of scattered 
light and/or diffuse emission when the nebula has a moder- 
ate to high absorption optical depth. For low optical depths 
re-absorption of diffuse emission in the nebula is negligible 
and the assumptions in CLOUDY work very well. However, for 
moderate to high optical depths the re-absorption of diffuse 
emission that is produced in the outer parts of the nebula 
and is radiated inwards becomes important. The assump- 
tions made in CLOUDY make it impossible to account for this 
energy source, nor for the amount of flux absorbed in these 
regions. Since the circumstellar envelope is optically thick 
in the AGB phase, our calculations always start shortly af- 
ter the transition from the AGB to the post-AGB phase is 
complete. 

For low temperature models the only source of diffuse 
light at optical and UV wavelengths is scattering of cen- 
tral star light by dust grains; for the highest temperature 
models bound-free emission also plays a role. Since for the 
highest temperatures the models have a low absorption op- 
tical depth, the bound-free emission causes only minor prob- 
lems and we can judge the quality of the models primarily 
by investigating the (wavelength-averaged) scattering opti- 
cal depth. An example of this approach will be shown in 
Section 0. 

In the calculations, the central star was assumed to emit 
as a blackbody. All the models were calculated for a distance 
of 1 kpc. 



3.1 Density profiles of the circumstellar shell 

The density profiles of the expanding shell are calculated 
for the homogeneous and spherical case. Using the mass loss 
prescriptions described above and assuming an outflow ve- 
locity they can be calculated for any moment in time t. 



stellar temperature as function of the time during the AGB, 
in order to compute the mass loss history usin g equ ation (0). 



p{r,t) = 



with 



Afwind(ioj) 



4TTr^(t,tcj)Vc^p{tcj) 



(8) 



r{t,tcj) = -R.(tcj) + (i - tej)Ucxp(icj) 

the distance from the centre of the star, 7?* the stellar ra- 
dius; tej stands for the time when a certain layer was ejected, 
iiexp is the expansion velocity of the wind at the moment of 
ejection. It is assumed to be constant thereafter and hy- 
drodynamical effects are neglected (see also Section 3.3). In 



particular, the post-AGB wind has a higher velocity than 
the AGB shell and will eventually overtake it. That part of 
the post-AGB wind which has done so contains little mass 
and is simply discarded. 



3.2 The AGB circumstellar shell 

The AGB shell is the principal contributor to the IRAS 
fluxes, and it is necessary to have a good description of the 



The Mcnv-Tf.ff relation discussed in Section 2.1.2 starts at 
temperatures roughly between 3000 K and 4000 K (the rela- 
tions we received from Blocker extend to lower temperatures 
than given in his paper). A description of the AGB evolution 
for lower temperatures is lacking. Unfortunately, the evolu- 
tion of Tcff and L on the AGB is not shown in the Blocker 
papers. We therefore simply extrapolate the Mcnv-T^s rela- 
tion logarithmically to lower temperatures. The 'AGB' star 
then evolves according to the extended relation. The M^nv- 
Tcff relation rises rather steeply at the low temperature end 
so that only a limited amount of extrapolation is necessary. 
In this way we find reasonable start temperatures for the 
AGB. Since our method implicitly assumes that the lumi- 
nosity remains constant during the AGB, the Bowen mass 
loss rates decrease with temperature. We investigate two 
different cases of AGB mass loss in the extrapolated part 
of the Mcnv-Tcff relation: the normal Bowen mass loss and 
a constant mass loss held at the value of the Bowen mass 
loss rate at the first point of the Mcnv-TcB relation as given 
by Blocker. These choices do not have implications for the 
post-AGB evolution. 



3.3 Expansion velocities 

The observed mean expansion velocity of AGB winds is 
15 kms~^ (Olofsson 1993), and this will be used as the typ- 
ical AGB outflow velocity. During the post-AGB phase the 
situation is different. Slijkhuis & Groenewegen (1992) as- 
sumed that the post-AGB wind has the same outflow veloc- 
ity (15 kms~^) as the AGB wind. The escape velocity (and 
hence also the outflow velocity) increases with temperature 
however. The increasing radiation pressure from the hotter 
star on the less dense post-AGB wind accelerates the dust 
and will decrease the densities in the post-AGB wind. The 
escape velocity for a 10 000 K star is already of the order of 
100 kms~^ to 150 kms~^. In addition, Szczerba & Marten 
(1993) found that dust was accelerated to 150 kms~^ dur- 
ing the post-AGB phase. We therefore assume an expansion 
velocity of 150 kms~^ in the post-AGB phase. 

One should realize that the scenario of a fast wind that 
follows a slow wind results in a collision between the two 
winds (as shown by e.g. Mellema 1993 and Frank et al. 1993). 
When inspecting the plots of Mellema (1993) we find that 
the effects of the colliding winds only start to become signifi- 
cant when the radiation driven wind (Mcpn) with velocities 
in excess of thousands of kilometers per second has devel- 
oped. The effect of colliding winds will be neglected in the 
further calculations in this paper. 



4 THE MODEL RUNS 

In this section we will investigate the spectral evolution of a 
post-AGB star by varying certain parameters that influence 
the mass loss rate, wind velocity and the evolutionary speed 
of the central star. We restrict ourselves to the 0.605 M© 
track and the emphasis will be on the infrared properties of 
the dust shell around the star. We calculated a total of seven 
runs. Every individual run consists of four different series of 
models. These individual models represent post-AGB shells 
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Table 3. The calculated runs for the 0.605 M© track, with an 
initial envelope mass of 2 Mq , and an AGB expansion velocity of 
15 kms~^. 



Run 


I'PAGB 


Pa 


Ph 


Constant 




(kms-1) 


(d) 


(d) 


mass loss 


1 


15 


100 


50 


no 


2 


15 


100 


50 


yes 


3 


150 


100 


50 


yes 


4 


150 


100 


50 


no 


5 


150 


125 


75 


no 


6 


1500 


100 


50 


no 


7 


150 


125 


75 


yes 



Table 4. The ages of the different models since the end of the 
AGB transition phase. The age of the AGB shell, in the case of 
the (Pa,Pb) = {100d,50d) models (= normal end), is 6465 yr in 
the case of a decreasing mass loss and 10 326 yr in the case of a 
constant mass loss. For the (Pa,Pb) = {125d,75d) models (=early 
end), these ages are 6741 yr and 10 603 yr respectively. 



with carbon-rich or oxygen-rich dust, both with and without 
dust formation in the post-AGB wind. We assumed the total 
mass of the AGB and post-AGB shell to be 2 M©. With the 
extrapolated Monv-Teff relation this implies a start temper- 
ature of the AGB of 2543 K. The different run parameters 
are outlined in Table w. 

The main differences between the runs are threefold. 
Firstly, as stated above, the post-AGB wind velocity is sub- 
ject to uncertainty. In order to assess the influence of this ve- 
locity we used three different values: 15 kms~^, 150 kms~^ 
and 1500 kms~^. Secondly, the differences between a de- 
creasing AGB mass loss rate and a constant mass loss rate 
are investigated. In the second case the constant mass loss 
is kept at the value of the Bowen mass loss rate at the 
first point of the Mcnv-Toff relation as given by Blocker. 
This value for the mass loss of Af = 1.99x 10"" M© yr"^ will 
then be adopted throughout the extrapolated part of the 
Monv-Teff relation, i.e. for temperatures between 2543 K 
and 3743 K. Approximately 85 % of the total shell mass 
is ejected in this phase. Thirdly, the definition of the end of 
the AGB is adjusted. As an alternative the pulsation periods 
that define the start of the transition wind and the start of 
the post-AGB phase are set to 125 d and 75 d respectively, 
which corresponds to a transition phase occurring between 
4723 K and 5418 K. For these runs the start of the transition 
period (defined by Pa) will be reached earlier, however the 
transition period itself will take much longer due to the fact 
that the Menv-reff relation is much steeper now in the tran- 
sition region, making the evolutionary speed much slower. 
This gives the paradoxical result that an earlier start of the 
transition phase gives rise to a later start of the post-AGB 
phase. Also the smaller mass loss rates result in a slower evo- 
lution between 5418 K to 6042 K (the start of the post-AGB 
phase in the (Pa, ft.) = (100d,50d) runs). After the star has 
reached Toff = 6042 K, the post-AGB mass loss rates are the 
same, and consequently the evolutionary speed in terms of 
kelvin per year is identical. 

The models are calculated at the temperatures listed 
in Table 0. The second and third entries in this table reflect 
the time that has passed since the end of the transition wind 
for the (Pa,Pb) = (100d,50d) and (Pa,Pb) = (125d,75d) case 
respectively. 

We will start with the results of run 4, which we will use 
as a reference since we consider the parameters for this run 
the most realistic. Subsequently we discuss the main differ- 
ences between this model run and the other runs. The main 
results of run 4 are visualized in Fig. M, where the tracks in 



Teff 


t(normal) 


t (early) 


Tcff 


t (normal) 


t (early) 


(K) 


(y) 


(yr) 


(K) 


(y) 


(yr) 


5 450 




553 


9 000 


904* 


5510 


5 500 




1308 


9 500 


926* 


5532 


5 550 




1950 


10 000 


948 


5553 


5 600 




2492 


11000 


991* 


5597 


5 700 




3328 


12 000 


1036 


5642 


5 800 




3904 


13 000 


1082* 


5688 


6 000 




4528 


14 000 


1131 


5736 


6 050 


13 


4619* 


15 000 


1180* 


5786 


6100 


89 


4695* 


16 000 


1232 


5838 


6150 


154 


4760* 


17 000 


1286* 


5892 


6 200 


211 


4817* 


18 000 


1342 


5947 


6 250 


262 


4867 


19 000 


1399* 


6004 


6 500 


457 


5063 


20 000 


1457 


6063 


7 000 


715 


5320 


25 000 


1768 


6374 


7 500 


823 


5429 


30 000 


2093* 


6699 


8 000 


860 


5466 


35 000 


2409 


7014 


8 500 


883t 


5488 









* These models are not shown in the graphite tracks in Figs. 
and m. 

* These models are not shown in the tracks in Fig. hol 



the IRAS colour-colour diagram are plotted, in Fig. where 
certain fluxes are plotted, and in Figs. H and H where a se- 
lection of spectral energy distributions is presented. In these 
flgures a distance of 1 kpc to the object is assumed. We will 
discuss the results for the oxygen-rich and carbon-rich dust 
separately. 

For this run we have made an estimate for the 
wavelength-averaged scattering optical depth (see also the 
discussion in Section H). For the models without post-AGB 
dust formation the results are shown in Fig. |j. The results 
for the models with post-AGB dust formation are not shown, 
but they are almost identical. It can be seen that the scat- 
tering optical depth for the silicate and graphite models are 
very similar, while on the other hand the absorption op- 
tical depth is much higher for the graphite models when 
compared to the silicate models (this can be judged by com- 
paring the V magnitudes for both models). This is caused 
by a combination of two effects. Firstly, graphite is a very 
efficient absorber at optical wavelengths and a relatively less 
efficient scatterer. The reverse is the case for silicates. This 
tends to level out the differences. Secondly, in the averaging 
process the optical depths are weighted by the output spec- 
trum. This spectrum peaks much more towards the red for 
the graphite models due to the higher internal extinction. 
The scattering optical depth is lower at longer wavelengths 
and this tends to level out the differences even further. The 
latter effect also explains why the scattering optical depth 
drops off so slowly, or even rises towards higher temper- 
atures. The peak of the energy distribution shifts towards 
the blue, where the scattering efficiency is much higher, thus 
countering the effects of the dilution of the circumstellar en- 
velope. 
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Oxygen rich 



Carbon rich 




60 /iin 

I Iiii 

10 20 30 10 20 30 

Figure 4. The V magnitude and the 12 fim, 25 fim and 60 ^m 
flux densities as function of temperature during the post-AGB 
evolution for the 0.605 Mq run 4 models. The soUd Hnes represent 
the models without dust formation in the post-AGB outflow. The 
dashed lines indicate the models with dust formation in the post- 
AGB wind. In the panels where no dashed hnes are visible, the 
values are indistinguishable. All curves start at T^ff = 6100 K, 
except the curves for graphite with post-AGB dust formation 
which start at 6200 K. 



Using these results we were able to obtain a worst case 
estimate for the amount of energy that could have been 
missed in the dust emission due to the fact that the scatter- 
ing processes were not properly treated in our models. The 
amounts for the models shown in Fig. H are typically 10 % to 
20 % of the total far-IR flux for the silicate models, and 5 % 
to 10 % for the graphite models. The worst case estimates 
for the models not shown in Fig. Ware all below 40 %. These 
estimates reflect on the accuracy of the absolute fluxes and 
magnitudes predicted by our models. However, they are ex- 
pected to cancel in a first order approximation when colours 
are calculated. Therefore we decided to omit the lowest tem- 
perature models whenever absolute fluxes are shown, but to 
include them when colours are shown. 

Scattering processes will not only influence the total 
amount of energy absorbed, but also the run of the dust 
temperature with radius. This can influence the shape of 
the spectrum and hence also the colours. It is impossible to 
make an estimate for the magnitude of this effect and we 
will assume it to be negligible. 
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4.1 Silicate dust 

We will start with a description of the silicate tracks in the 
IRAS colour-colour diagram (Fig. |). At first the [12]-[25] 
and [25]-[60] colours both increase. At some point the [25]- 
[60] colour remains constant, while [12]-[25] still increases. 
This is followed by a decrease in both the [12]-[25] and [25]- 
[60] colours. Hence the track makes a clockwise loop in the 
colour-colour diagram. 

The first part of the track represents the cooling of the 
dust shell due to its expansion. From Table ^ one sees the 
inner radius of the shell has increased by a factor of 4 be- 
tween Toff = 6200 K and T^a = 8000 K, which leads to 
a cooling and a decrease of the IRAS fluxes (Fig. M). The 
following evolution in the IRAS colour-colour diagram is 
rather counter-intuitive. Normally one would expect that 
the shell would make the familiar counter-clockwise loop in 
the colour-colour diagram as found by Loup (1991), Volk & 
Kwok (1989) and Slijkhuis & Groenewegen (1992). In such 
a loop, the shell continues to cool, until the photospheric ra- 
diation begins to dominate the emission, first at the 12 /im 
band. At this time, the 25 /im flux and the [12]-[25] colour 
decrease strongly due to the expansion and cooling of the 
shell. The [25]-[60] colour will continue to increase slowly. 
Later, as the star starts to dominate at 25 /xm, the [12]-[25] 
colour will remain constant and the [25]-[60] colour will start 
to decrease. When the star is the dominant contributor in 
all three IRAS bands, the loop will end in the Rayleigh- 
Jeans point. However, the above authors assumed a con- 
stant temperature of the central star, while the effects of an 
evolving star on the IRAS colours can not be neglected; in 
the 0.605 M© track the star rapidly evolves toward higher 
effective temperatures. For example, the kinematic age of 
the shell, and thus its outer (or inner) radius, has increased 
by only 25 % between T^s = 8000 K and T^s = 12 000 K. 
One can regard the circumstellar shell as essentially sta- 
tionary around the evolving star. An evolving star embed- 
ded in a stationary shell results in a heating of the dust 
and an increase of the IRAS fluxes (Fig. W). The effect of 
re-heating of the circumstellar envelope was already found 
by Marten, Szczerba & Blocker (1993) in their calculations. 
The re-heating continues even beyond the last data point, 
hence in the colour-colour diagram the track keeps evolving 
to higher colour temperatures in both colours. Around the 
turning point the 12 /im flux is partly due to the photo- 
sphere. The photospheric flux decreases for higher temper- 
atures so that the 12 pm flux reacts later to the rising dust 
flux than the 25 fj,in and 60 pm flux. Therefore the [12]-[25] 
colour starts decreasing later than the [25]-[60] colour. The 
result is a clockwise loop in the colour-colour diagram. 

One should note the behaviour of the IRAS flux den- 
sities in Fig. W. The flux densities reach a minimum during 
the slow evolution of the star. The minimum is followed 
by gradual increase in the total infrared output due to the 
heating of the shell resulting from the increasing absorbing 
efficiency of the dust (see also the discussion in Section 4.2). 



A model star like this (i.e. with silicate dust) would have 
a larger chance of being detected in an infrared survey like 
IRAS when it has evolved to higher temperatures. This fact, 
combined with the predicted distribution over spectral type 
implies that hot post-AGB stars, or young PN should be 
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- ^ - graphite 
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[12]-[25] 

Figure 3. IRAS colour-colour diagram — run 4. The arrows indicate the evolution of the IRAS colour with time. The definitions for the 
IRAS colours are as presented in Oudmaijer et al. (1992). With these definitions the Rayleigh- Jeans point, where objects hotter than 
approximately 2000 K are found, is located at (0,0). The boxes as defined by van der Veen & Habing (1988) are drawn in. The straight 
line represents the blackbody line. Selected temperatures are indicated. The first data point for all tracks is the T^ff = 6050 K model. 



expected to be more abundant in samples of evolved stars 
with infrared excess. 



track will move asymptotically to the track without post- 
AGB dust formation. 



4.1.1 Dust formation in the post-AGB wind 

The addition of dust formation in the post-AGB wind 
changes the spectral energy distribution. From Fig. Ml it is 
visible that the 25 /im and 60 fim flux densities have the 
same values as in the AGB-only dust formation case. The 
12 pm band is strongly affected by the addition of the hot 
dust. This is due to the presence of the 10 ^m silicate emis- 
sion feature which causes the emissivity of the hot silicate 
dust to peak strongly at 10 fim. The larger 12 fim flux den- 
sity is immediately reflected in the colour-colour diagram. 
It puts the starting point at a slightly higher colour tem- 
perature than before. As the feature increases in strength, 
the track bends toward higher [12]-[25] colour temperatures. 
Then, when the star is evolving more rapidly than the shell 
expands, the cool AGB shell is heating up, and the track 
moves down. At the same time, the [12]-[25] colour temper- 
ature decreases, reflecting a weakening of the silicate fea- 
ture relative to the 25 fim flux. This effect is caused by the 
circumstances under which dust in the post-AGB wind is 
formed. As the star becomes hotter, also the distance in- 
creases at which the equilibrium temperature of the dust 
goes below 1000 K. Therefore the dust formation will take 
place at larger distances from the star, in a diluted local 
radiation field where the dust density will be lower. Conse- 
quently, the silicate feature weakens and the track moves to 
lower [12]-[25] colour temperatures. Eventually the contri- 
bution of the 10 fim feature will become negligible and the 



4.1.2 Different post-AGB wind velocities 

In Fig. Vn the effect of a different post-AGB wind velocity is 
shown for silicate dust in the upper panel. The models have 
been calculated for velocities of 15 kms~^, 150 kms~^ and 
1500 kms"\ The 1500 kms"^ tracks resemble the AGB- 
only dust formation models, because the amount of dust in 
the post-AGB wind is not large enough to yield an observ- 
able effect. The low outflow velocity model shows all the 
effects outlined for the 150 kms~^ models even stronger be- 
cause the amount of hot dust has increased. 



4.2 Carbon-rich dust 

Let us now turn to the models with carbon- rich dust. This 
type of dust is a more efficient absorber than silicate dust as 
can be seen in Fig. H. The visual magnitude is larger than 
for the oxygen rich models, reflecting a larger circumstellar 
extinction in the Johnson V band. One can also see that the 
infrared energy output is larger. 

The beginning of the track of the carbon-dust mod- 
els in the colour-colour diagram (Fig. m is located at a 
higher colour temperature than for the oxygen-rich models 
because more energy is absorbed by the dust. At first, the 
track moves to lower temperatures, and as the star begins 
to evolve more rapidly, the curve makes a slight backward 
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Figure 7. In these //?j4 5 colour-colour diagrams, the effect of different velocities of the post-AGB wind with dust formation is presented. 
The 15 kms-^, 150 kins-^ and 1500 kms~^ models are calculated in runs 1, 4 and 6. The upper diagram shows the silicate models, the 
lower diagram the graphite models. The graphite 15 kins^^ track starts at 6200 K, all others at 6050 K. 
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Figure 5. The spectral energy distribution at several stages dur- 
ing the post-AGB evolution. The upper panel shows the SED for 
silicate dust without post-AGB dust formation, the lower panel 
with dust formation, both for run 4. The flux densities are cali- 
brated for an assumed distance of 1 kpc. The temperatures are 
indicated. Note the presence of the silicate emission feature at 
10 /^m in the post-AGB dust models. 



loop in the diagram due to the heating of the shell. Then, 
contrary to the silicate model, the shell cools again. 

This somewhat unexpected result can be explained by 
the properties of the absorbing material. It is instructive to 
investigate the effective cross section Qo of the two grain 
species as a function of the effective temperature of the cen- 
tral star spectrum (i.e. a blackbody). Qo is defined as 



Qo 






(9) 



With Bu the blackbody intensity distribution and a„ 
the absorption cross section of the grains. The formula is 
chosen in such a way that the total energy absorbed by the 
grains is proportional to L x Qo- Note that L is constant 
during the evolution. The absorption coefflcients have been 
chosen such that both grain species have the same dust-to- 
gas mass ratio. The resulting curves are shown in Fig. 

For temperatures between approximately 1000 K and 
50 000 K silicates absorb less energy than graphite. It can 
be a factor of 10 lower in the temperature range for cool 
post-AGB stars. This is caused by the fact that silicates are 
inefficient absorbers at optical wavelengths. From these re- 
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Figure 6. As the previous figure, but now for graphite with and 
without post-AGB dust formation (lower panel, respectively up- 
per panel). Note the hot dust that emits strongly in the near- 
infrared, but hardly contributes to the 12 /im flux density. 
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Figure 8. The effective cross section for the grain species used 
in our study as a function of the effective temperature of a black- 
body. The total luminosity of the central source remains constant 
and the dust-to-gas mass ratio for both grain species is equal. 
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suits we can expect that for temperatures between 10 000 K 
and 50 000 K the total IR emission of silicates will rise much 
more drastically than for graphite. For temperatures above 
50 000 K the peak of the energy distribution shifts into the 
EUV. At these wavelengths silicates are more efficient ab- 
sorbers and we can see that for these temperatures the total 
amount of energy absorbed by the silicates is higher than 
for graphite. 

Thus, in the last data points of the track in the IRAS 
colour-colour diagram, the absorption efficiency of graphite 
does increase less strongly with rising effective temperature. 
In this case, the interplay between the evolving star (now 
heating up the shell less rapidly) and the expanding shell 
(giving rise to cooling) has as a result that the expanding 
shell becomes the more dominant factor and consequently 
the shell cools. In Fig. W these two effects are illustrated by 
the decreasing IRAS fluxes. 



4-2.1 Dust formation in the post-AGB wind 

In contrast to oxygen-rich dust, graphite has no feature in 
the IRAS 12 pm pass band. Also the shape for the emissivity 
law of graphite is such that the hot dust peaks at near- 
infrared wavelengths, contributing negligibly to the 12 fira 
band (Figs. H and pi). Therefore one does not expect a large 
effect of the hot dust on the IRAS 12 fim flux. However, the 
12 /xm, 25 fim and 60 ^m flux densities of these models are 
lower than for the corresponding AGB-only dust formation 
models (Fig. W. This stems from the fact that the hot dust 
absorbs a large fraction of the stellar energy, yielding a cooler 
AGB shell. 



4-2.2 Different post-AGB wind velocities 

The effect of the post-AGB outflow velocity on the 
graphite models is presented in the lower panel of Fig. \A 
The 1500 kms~^ track behaves almost the same as the 
150 kms~^ track and the track without post-AGB dust. 
The 15 kms^^ track however is different, first it moves up- 
ward, then bends downward and finally makes the counter- 
clockwise loop. The initial increase in the [25]-[60] colour is 
caused by the fact that the large amount of hot dust now 
shields the AGB shell more efficiently, making the AGB shell 
even cooler than previously. As the star becomes hotter, the 
dust condensation radius moves rapidly outward, making 
shielding less efficient. The AGB shell now obtains more en- 
ergy from the central star and heats up, resulting in higher 
[25]-[60] temperatures. The tracks evolve asymptotically to- 
wards each other and finally the counter-clockwise loop sets 
in, as was already explained. 



4.3 Constant mass loss 

In Fig. the silicate tracks for run 4 (decreasing AGB mass 
loss) and run 3 (constant mass loss) in the IRAS colour- 
colour diagram are shown. Both the [12]-[25] and the [25]- 
[60] colour temperature of the constant mass loss track are 
slightly higher than for the decreasing mass loss track. 

The constant mass loss creates relatively more cool dust 
further from the star. This cooler dust emits less radiation 
and therefore the hotter dust at the inner parts of the AGB 



shell (where the mass loss rates for both tracks are the same) 
will dominate the spectrum more. This gives the result that 
the integrated spectrum of all dust appears hotter. This ef- 
fect is strongest at 60 fim, less at 25 fim and absent at 
12 ^m. As a whole, the effect of different AGB mass loss 
rates is small. The effect on the graphite tracks is basically 
the same and is not shown. 



4.4 Different end of the AGB 

For all the models that have been discussed so far, it was 
assumed that the transition mass loss starts at Pa ~ 100 d, 
and that the Reimers mass loss starts at Pb = 50 d. How- 
ever, the evolutionary behaviour of the calculated models 
strongly depends on the exact moment of 'superwind' ces- 
sation (Szczerba 1993). In this section, we will present an 
investigation into the effect of changing this unknown pa- 
rameter. 

As stated before, the evolution of the central star in 
the (Pa,Pb) = (100d,50d) and the (Pa,Pb) = (125d,75d) 
models is identical after Tcs ~ 6042 K. However, the den- 
sity structure of the circumstellar shell is entirely different 
when the star has reached this temperature. Whereas the 
(Pa,Pb) = (100d,50d) model then barely has a detached 
shell, the (Pa,Pb) = (125d,75d) model already has a shell 
with a kinematic age in excess of 4500 yr, which has much 
cooler dust. 

We will now look at the results for this run in the colour- 
colour diagram. It appears that the alternative end of the 
AGB entirely changes the path in the colour-colour diagram 
(Fig. fU]) . A selection of spectral energy distributions is pre- 
sented in Fig. [01. 

From the start, the silicate dust track moves strongly to 
the upper left, making a turn to the lower right, and finally 
a bend to the left. Since the star evolves slowly (it takes 
approximately 750 yr to go from 5450 K to 5500 K), the 
circumstellar shell cools rapidly, and only the stellar pho- 
tosphere is visible at 12 pm (Fig. hlh. This is reflected in 
an increasing [12]-[25] colour temperature: the photospheric 
12 ^m flux remains approximately the same, while the 25 ^m 
flux and the [25]-[60] temperature decrease due to the cool- 
ing of the shell. Subsequently the star starts to evolve more 
rapidly than the shell expands, and the [12]-[25] colour tem- 
perature decreases. This is due to the fact that the shell 
now heats up, increasing the 25 /im flux and decreasing the 
[25]-[60] colour, contrary to the 12 /im band which is still 
dominated by the decreasing photospheric emission of the 
evolving star. The shell continues to heat up, and the dust 
emission at 12 /xm overtakes the decreasing photospheric 
flux density. The curve bends to higher [12]-[25] tempera- 
tures. The contribution of nebular bound-free emission at 
12 fim for the hottest models increases the [12]-[25] colour 
temperature even more, and the track makes the, now fa- 
miliar, clockwise loop. 

All effects outlined above appear less strong for the 
carbon-rich model. It initially starts on a cooling curve, and 
later bends to the clockwise loop. This is explained by the 
fact that the dust emission from the carbon grains still is 
present at 12 pm in the first points, so that the shell cools 
in all colours. For the remainder of the evolution there al- 
ways remains some dust emission in the 12 /xm band, making 
the loop to the left less pronounced. 
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Figure 9. IRAS colour-colour diagram in which the difference between the decreasing and constant mass loss (runs 4 and 3) arc shown 
for the silicate tracks. The first data point for both tracks is the T^g = 6050 K model. 
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Figure 10. IRAS colour-colour diagram for the models without post-AGB dust formation in run 5, the case with an earlier start of the 
transition phase than in run 4. The first data point for both tracks is the T^ff = 5450 K model. 



© 0000 RAS, MNRAS 000, 000-000 



The spectral evolution of post-AGB stars 15 





3 








- 








/■',-' — ~.v. 




- 




2 




■■/?^^. 




- 








.' / / ■■• \\ 




- 




1 


RX . 


/;/ \ 




- 




n 




/' .■ 5450 K '\'^ 

' / 6500 K Y^ 




-■- 








,' ■■ 12000 K ^ 

20000K 

35000 K 


N^x,^._ 


- 




-1 


- '"■'''^'' 




'\ 


- 



1 2 3 

lg(?L/|im) 



c' 1 







/ /'f 


*^^ 










/ ■'." ■' 


\ 




. 


- 




■V' ■ 

//'.■■' 
//'.■■ 


\ 




- 


-^ - — 


/ 

/ / 

^\_y ill 


/■■■ 

'',■■ 


5450 K 

6500K 

12000 K 

20000K 

35000K 


\ 


^:] 


^ 1 , , 


1 




1 




\ 



12 3 
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without post-AGB dust formation. The flux densities are cali- 
brated for an assumed distance of 1 kpc. The temperatures are 
indicated. 



gram. We have only investigated this for the 3 M0 track 
and it is not yet clear if this observation is valid in a more 
general context. In the past years the /i? A 5 colour-colour di- 
agram has proven to be a very useful tool for studying post- 
AGB evolution. Our knowledge has increased considerably 
since its introduction. However, when we try to understand 
the details of this evolution better, the information from the 
IRAS colour-colour diagram becomes more and more con- 
fusing. The K-[12] vs. [12]-[25] diagram may be a valuable 
additional tool for the study of post-AGB evolution. 

We will only discuss the evolution of the -K'-[12] colour in 
Fig. h2l since the evolution of the [12]-[25] colour already has 
been discussed earlier. At first both the fiux in the K band 
and the 12 fim band are decreasing. However the 12 fim flux 
decreases more rapidly and therefore -K'-[12] evolves towards 
hotter colours. This can be understood if we realize that the 
central star heats up relatively slowly while the circumstellar 
shell expands relatively rapidly. At a temperature of around 
8000 K the central star evolution will speed up considerably, 
reversing the preceding argument. The flux in the K band 
continues to drop, while the flux in the 12 /im band increases 
now, resulting in ever cooler K-[12] colours. This evolution 
continues until the central star starts to ionize a considerable 
part of the circumstellar shell and the K band flux will start 
to rise again due to nebular bound-free emission. This rise 
is so rapid that it reverses the evolution of the -K'-[12] colour. 

The tracks for silicate and graphite grains are quali- 
tatively similar, but are offset with respect to each other. 
This is mainly due to the different absorption efficiency of 
silicate and graphite grains (see Section 4.2). The difference 



between the silicate track with and without post-AGB dust 
formation can be understood solely from the presence or ab- 
sence of the 10 fim emission feature. The difference between 
the graphite tracks with and without post-AGB dust forma- 
tion stems from the fact that the hot dust in the post-AGB 
section of the wind mainly radiates around 3 /im and thus 
contributes to the K band flux. 



4.5 A near-IR colour-colour diagram 

In this section we present the evolution of a post-AGB star in 
a colour-colour diagram which uses colours at shorter wave- 
lengths: the K-[12] vs. [12]-[25] diagram. This diagram gives 
a different view of the evolution since in all models the John- 
son K band will be dominated either by the stellar contin- 
uum or the bound-free emission of the ionized part of the 
nebula. Only in the models containing hot graphite grains, 
part of the flux in this band will originate from the grains. 
So, contrary to the IRAS colour-colour diagram which we 
presented earlier, this diagram contains information both 
on the central star and the dust. The results for run 4 are 
shown in Fig. Fl2, giving both the tracks using silicate and 
graphite grains and both the tracks assuming AGB-only and 
post-AGB dust formation. 

The flrst thing we notice is that all four tracks, at least 
in a qualitative sense, look quite similar. This suggests that 
this colour-colour diagram is less sensitive to particulars of 
the grain emission and thus that the information on the 
evolution of the central star and the nebula is less 'contam- 
inated' when compared with the IRAS colour-colour dia- 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have presented a new model to calculate 
the spectral evolution of a hydrogen burning post-AGB star. 
The main new ingredient of this model is the possibility of 
extracting timescales of the post-AGB central star evolution 
from the most recent evolutionary calculations. Hence, con- 
trary to previous studies, it is possible now to investigate 
other AGB and/or post-AGB mass loss rates and different 
prescriptions for the start of the post-AGB phase. The use 
of a photo-ionization code in which a dust code is built in, 
gives us the opportunity to study the evolution of the in- 
frared emission of the circumstellar material. 

We have performed a parameter study on a typical post- 
AGB star with a core mass of 0.605 M© taken from Blocker 
(1995b). By varying the parameters that govern the mass 
loss and the evolutionary timescales of the post-AGB phase 
by a reasonable amount we flnd that: 

(1) The influence of the evolving star on the infrared 
colour evolution can not be neglected. In particular, the 
phase wherein the circumstellar shell can be considered as a 
stationary shell around a star that rapidly increases in tem- 
perature results in clockwise loops in the /i? AS colour-colour 
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Figure 12. The ncar-IR colour-colour diagram at several stages during the post-AGB evolution for run 4. For reference a blackbody 
curve is also included in the diagram. 



diagram. This was not observed in many previous studies in 
which the temperature of the central star was taken to be 
constant. Instead, in these studies the evolutionary path fol- 
lowed a counter-clockwise loop in the colour-colour diagram 
because of the inevitable cooling of the shell. 
Only Szczerba & Marten (1993) who used a Blocker track 
obtained a roughly similar result with their dust radia- 
tive transfer code. Volk (1992) found a slight deviation 
from the counter-clockwise loop in the colour-colour dia- 
gram. He used the coarse grid of evolutionary timescales 
from the 0.644 M© Schonberner (1983) track. The output of 
the photo-ionization code CLOUDY was used as input for a 
dust radiative transfer model. It is not clear to what extent 
that choice influenced the path followed in the /i?y4 5 colour- 
colour diagram. 

The crucial influence of the evolution of the central star 
warrants further parameter studies with other stellar mod- 
els, where the increase in temperature as a function of time 
will be different. It is expected that the heating of the shell 
will be more important for more rapidly evolving (i.e. larger 
core mass) stars. 

(2) The tracks that are followed in the /i?ylS colour-colour 
diagram are very sensitive to the adopted dust opacity law 
and solid state features. Two main differences between the 
silicate and graphite models that were studied are evident. 
Firstly the IRAS colours of the silicate models are very sen- 
sitive to newly synthesized dust in the post-AGB wind be- 
cause of the silicate 10 pm feature that contributes signifi- 
cantly to the 12 /jm flux density. In contrast, hot graphite 
dust emits mainly shortward of the 12 /xm pass band and 
thus post-AGB dust formation has less influence in this case. 
Secondly, the dependency of the absorption efficiency on the 



central star temperature is much stronger for silicates than 
for graphite in the temperature regime studied here. There- 
fore silicates react much stronger to the heating of the cen- 
tral star and thus give rise to much larger loops in the IRAS 
colour-colour diagram. 

Our knowledge of dust opacities, and certainly of solid state 
features in the mid- and far-infrared, will improve when the 
results of the 750 mission have been digested (Waters et 
al. 1996). For example, the well-known 21 fim and 30 ^m 
features that have been observed in the infrared spectrum 
of carbon-rich post-AGB stars and planetary nebulae (e.g. 
Omont et al. 1995) have not been taken into account in this 
study. In addition, the wavelength coverage up to 200 fim 
will be of great help to determine the wavelength depen- 
dence of the dust opacities towards long wavelengths. 

(3) A third decisive factor that governs the evolution of 
the IRAS colours is the definition of the end of the AGB. The 
sooner an object enters the transition phase, i.e. the sooner 
the heavy AGB mass loss ceases, the longer the evolution to 
higher effective temperatures will last. This results in cool 
circumstellar dust shells, and consequently these models are 
on a location in the IRAS colour-colour diagram where not 
many post-AGB stars were expected previously. In this re- 
spect it is noteworthy to refer to Fig. nO where oxygen-rich 
post-AGB stars are predicted to be present in the upper 
part of region VIII. A re-investigation of the sources in that 
region would be useful in testing whether this scenario is 
realistic. 

(4) Changing the mass loss prescription in the coolest 
part of the AGB evolution that we considered has little in- 
fluence on the IRAS colours of the models. 

In general we find that the variation of the parameters 
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mentioned above, which are stiU not very well determined, 
result in a variety of different paths in the IRAS colour- 
colour diagram. This is certainly part of the explanation 
why planetary nebulae do not occupy a well-structured re- 
gion in the fflAS colour-colour diagram (cf. Volk 1992). As a 
by-product of this investigation we find that the same loca- 
tion in the IRAS colour-colour diagram can be occupied by 
objects with an entirely different evolutionary past. Appar- 
ently, the location in the IRAS colour-colour diagram can 
not a priori give a unique determination of the evolutionary 
status of an object. 

As an alternative to the IRAS colour-colour diagram, 
the -K"-[12] vs. [12]-[25] colour diagram is presented. The 
tracks in this diagram seem less affected by particulars of 
the grain emission. Hence this diagram might prove to be a 
valuable additional tool for studying post-AGB evolution. 

The feedback of observational work, after a sufficient 
number of parameter studies, will be of help to assess the 
selection effects in our post-AGB sample selection criteria, 
and will give constraints on the assumptions that now have 
to be made on the central star and nebular evolution. 
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